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1. INTRODUCTION 

This is the first semi-annual report on a study of Gimballess 
Inertial Navigation Systems. All of the tasks' yhich were originally 
proposed for inclusion in this study -are being actively pursued. 

They are: (l) a comparative study of several methods of measuring 

angular velocity without the use of gyroscopes, (2) a study of methods 
of error damping for flights that are not earth-hound, (3) the 
construction of a breadboard model of a rotating-accelerometer device 
for measuring angular velocity, (k) the determination of the computer 
requirements for performing the computations necessary for gimballess 
inertial navigation systems, (5) the simulation of one or more mechaniza- 
tions of gimballess systems and (6) an extensive error study of gimballess 
systems . 

A su m mary of the work of the first six months is given in the 
next section of this report and details are given in the appendices . 

The technical staff has consisted of two Research Fellows, each 
at 75$ of full time, and one Associate Professor at 20$ of full time. 

In addition, there has been part time programming assistance. It is 
anticipated that by the end of the first year of this study, both of 
the Research Fellows will have completed all of their work for the Ph.D. 
degree and that further work on gimballess inertial navigation systems 
would be carried out by other advanced graduate students . 
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2. SUMMARY OP -WORK 

2.1 Measuring Angular Velocity -without Using Gyroscopes 

Most Inertial navigation systems use accelerometers to sense 
linear acceleration and gyroscopes to sense angular velocity. (Che use 
o^f gyroscopes may be avoided, however, hy using six or more accelerometers 
that are fixed to the vehicle or hy mounting an accelerometer on each 

I / 

of two perpendicular rotating rings. 

Many configurations of fixed accelerometers will allow the 
determination of hoth linear acceleration and angular acceleration. 

Several of these configurations are discussed in some detail in 
Appendixes A and B . The minimum number of fixed accelerometers for an 
all -accelerometer inertial navigation system is six. t With six accelero- 
meters there is an ambiguity in the sign of the angular Velocity. This 
ambiguity may he. avoided with eight or nine accelerometers. Also, the 
use of eight, or more accelerometers provides greater freedom in the 
placement of the accelerometers. It is shown in Appendix B that the- 
linear and angular acceleration can he readily determined even though' the 
accelerometers are not placed symmetrically about the vehicle’s center of 
mass and even though the center of mass may move due to fuel consumption. 

Section 2.3 discusses the use of accelerometers mounted on rotating 

rings . 

2.2 Error Damping 

1 2 

Bodner and Seleznev and Krishnan have studied the mechanization 
of the navigation equations for gimballess navigation systems and have 
■ found that mechanization in an inertial frame of reference is more 
suitable than mechanization in a vehicular frame of reference. 
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However, a simple mechanization of the equations in an inertial frame 
gives a j characteristic equation of the form 

(s 2 + (b 2 ) (s 2 - 2 ^) = 0 

corresponding to a transient error that contains a hounded sinusoidal 
term and an unbounded hyperbolic cosine term. Appendix E treats three 
methods of achieving a stable mechanization of the navigation equations . 
They are (l) damping by means of reference trajectory information, 

(2) damping by means of external velocity information and (3) elimination 
of diverging errors by means of _ altimeter information . 

For many flights the reference (desired) trajectory of the vehicle 
is accurately specified and the actual trajectory does not deviate far 
from this reference trajectory. With the. reference trajectory informatior 
available on board the vehicle, damping can be obtained without' the use 
of any auxiliary sensors; the accelerometers required for an undamped 
system suffice . The accuracy of a system with reference trajectory 
damping is dependent upon the closeness of the actual and reference 
trajectories . 

Doppler radar has been used to obtain velocity information for 
use in airborne terrestrial navigation systems . For space navigation , 

3 

Doppler radar is less 'satisfactory. Franklin and Birx have reported 
encouraging results, from feasibility and accuracy . studies of optical 
Doppler. Doppler -damped navigation system equations for space flight 
are presented in Appendix E. 



The distance from the vehicle to a celestial body may be found 
by a variety of distance measuring schemes. The measured value of this* 
distance may be used to eliminate the diverging error that occurs in 
the undamped navigation system. However, it will not eliminate the 
sinusoidal oscillation. 

2.3 Rotating Accelerometer Experimental Studies 
2 

Kfishnan has’ shown that two or three linear accelerometers 
mounted on mutually perpendicular rotating rings can be used to determine 
both linear acceleration and angular velocity. Ajco.t&ting disk, mounted 
on a dividing head and carrying an accelerometer, is under construction. 
The experimental study -to determine the feasibility of a rotating aecelero 
meter sensing system has not yet been initiated due to delays in 
obtaining transformers for the accelerometer "sp-loop-closin^ 1 amplifier 
and the difficulty in obtaining satisfactory operation of the loop- 
closing amplifier . It is expected that the use of very carefully matched 
diodes will remove the remaining difficulty in the loop-closing amplifier 
dnd that experimental studies will commence -during May. 

2.k The Oomputer Simulation 

The navigation system is to be simulated - on a digital computer . 

This simulation will permit a more adequate analysis of the system 
than is possible by purely analytic means. It will also permit a. study 
of the navigation system errors and will provide -the necessary information 
for determining the requirements for an on-board computer. 
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In the simulation of the navigation system, the inputs are 
an angular velocity pattern and a linear acceleration pattern. 

With these inputs, the simulated system.- and the actual system - 
will provide the vehicle's position coordinates, which due to errors 
in the sensors, errors in the calculation of the gravity compensation, 
errors in the direction cosine calculations and errors in the numerical 
integrations, will deviate from the true value. 

Figure 1 shows the block diagram for the simulation of the 

navigation system on a digital computer. The linear acceleration pattern 

a^.(t) enters in the upper left hand corner of the figure. It is 

integrated twice, combined with initial conditions and the gravity .terms 

-i 

and fed to Box 1 where the exact direction cosine matrix [D(T. )]„ is 
formed. This direction cosine matrix relates the inertial coordinate 
system and the vehicular coordinate system. The inputs to Box 1 are 
referred to the inertial coordinate system and the outputs, are referred 
to the .vehicular coordinate system. Box 2 combines, the linear accelera- 
tion Ay(t), the angular velocity a^(t) and the acceleration r^(t) of 
the origin of the vehicular coordinate system [V] with respect to the' 
origin of the coordinate system [m] about which the accelerometers are 
centered. The outputs of Box 2 are the accelerations that are read .by 
the accelerometers. Box 3 simulates the accelerometers. Its inputs 
are the true accelerations and noise and its outputs are the accelerometer 

readings.. ' The calculated values of- the angular velocity and of 
_ c 

linear acceleration are formed in Box h. Since -the accelerometer 
c 
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readings are noisy a linear filter and estimator is included (Box 5). 

Box 6 forms the calculated direction cosine matrix [D(T. )3 The 

i c 

inputs to Box 6 are referred to the vehicular coordinate system 

and the output iL. (T. ) is referred to the inertial coordinate, system. 

This linear acceleration, (T^), is- combined with the calculated 

c 

gravity term gj. and the resultant is integrated twice to give the 
c 

vehicle's calculated position'll. (T. ) . A linear predictor (Bdx 12) 

I c 1 ' 

gives the estimated position at the next sampling instant, T^ +1 ,and 
the gravity computer (Box 8) gives the corresponding acceleration Of 
gravity. Boxes 13 and 1^ give the calculated and exact Euler Angles 
respectively . 

Portions of the ■ simulation have been successfully tested. 

The exact form of some boxes, such as the linear filter and, estimator, 
have not been determined, lior has a closed loop test .been attempted; 
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APPENDIX A 

ALL ACCELEROMETER TECHNIQUES FOR MEASURING 
LINEAR ACCELERATION AND ANGULAR ACCELERATION 

A. R. Schuler 

Currently, most inertial navigation systems use linear 
accelerometers .to sense linear accelerations and gyroscopes' to sense 
angplar velocity or angular position. It is possible, however, to 
determine both linear acceleration and angular acceleration (or 
anguilar velocity) without the use of gyroscopes. This report 
presents four configurations of linear accelerometers which permit 
the, determination of both linear acceleration and angular acceleration 
(or angular velocity). Although there are many possible configurations, 
only those that seem to be of major importance are treated here. It 
is to be noted that the possible configurations differ in the number 
of accelerometers required and the mathematical form of the outputs . 

In' this. section, it is assumed that the accelerometers are placed on 
a coordinate system that has as 'its origin the center of gravity of 
the vehicle. 

Two coordinate systems are used: 

a) the inertial system [i], and 

b) the vehicular coordinate system [v] which has its origin 
at the center of mass of the vehicle and its axes V^, Vg 
and Vj along the principal axes of the vehicle. 



9 - 


Figure A-l shows the coordinate system [v] which, in general, 
is rotating and .accelerating with respect to the inertial coordinate 
system [i]. The angular velocity of [v] is' designated as: 

3 = \^ + \^ + % 3 \ {A ^ 

It has been shown -1- that the - inertial acceleration of', an arbitrary 
point P accelerating with respect" to a moving reference: frame is given- 
by (a dot represents differentiation wijfch respect 'to. time): 

Aj = Rj + 8 x + 2 S - x + B x (d s r^) (A-2) 

where : 

Ej is the vector from the origin of the inertial frame to 
the origin of the vehicular frame - , 
r^ is the vector from the origin of the vehicular frame 
to the point P. 

If the point P is fixed in the vehicular system then 
r„ = x^. = 0 and the acceleration of the point is given by: 

Aj. = Bj+Qxr^+J2x(flx r ) (A-3) 

Rj may be written as a sum of components along the vehicular axes, ,i.e.: 

*1 = \ 


(A-4) 




.X 


!0xen 

R I ^ \ + '®V 2 + ®V 3 \ 

k v 1 \ + \ 2 -fy + \ 3 = \ ^ A " 5 ^ 

Hote that E^. is an inertial acceleration whose components 

have been resolved along the • vehicular system axes . We introduced 

the notation R = A.. (i = 1, 2, 3) . 
i i 

The acceleration in inertial space (A_j.) will now he found 
for the six points, a, b, c, d, e, and f, fixed on the vehicular 
axis as shown in Fig. A-2. 

To find the acceleration of point e with respect to inertial 
space, let 

= *e \ (A ' 6) 

Using Eqs. (A-l), (A-5), and (A-6) in Eq.. (A-3), 

®I^pt e = A V 1 *V + -V + \ 3 \ + ^ *V + *V g •V' + A e \ 

+ \ \ + \ X + % V X [ <\ \ + \ X + ^) X *eV 

= V \ + ^e % + *e \ 5 V 3 ] 

+ Jv \ - *e \ + 4 fc \ % ( A ‘ 7 > 

+ 1 V ^ ^ * e ^ 
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To obtain the inertial acceleration at the other fire points, 
let r^ take on the aesirea cuumrna.be vacua . 

For point f let r ? = - Jj f k^ 

^t) pt f = ^ “ % 2 & f " fl v 3 ■* 

+ 3y \ + \ *t - \ \ V U " 8) 

+ \ [ \ 3 + X *t + ^ V 

For point a let r^. = 

pt a = \ X x " °f 2 £ a " ^ ^a 3 

+ jy [A V 2 + % *a + \ \ *a ] ' A “ 9) 

+ ^ [J V 3 - \ ,^a + \ % *a 3 

For point h, let r^ = - i^ 

(A i } pt b = Sr \ + ^ 2 % + V 

+ % l X ' X ‘ (A " 10 

+ \ + \ *b “ \ 
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For point e let r = A 
y c V 

pt c = 1 Vj^ ' ■ 8 o' + ' J V 1 % 2 ^ 

+ V I\ " ^ *c “ ^ *e ] 

+ \ [ % + \ *c + \ \ X 

For point d let r y = -jJ d 3^ 

pt a = *v X 1 + ®v 3 % " X ^2 ^ 

+ \ + \ % + 

+ *\ 3 ' n v 1 % ” ^ ® f 3 %- 

Also, if r =0 

v 

origin = + 'V \ g ^ 


Consider first the configuration of linear accelerometers 


(A -11) 


(A-12) 


(A-13) 


shown in Fig, A-3 (configuration A). The six accelerometers are 
oriented about the vehicle *s center of gravity and are mounted on the 
principal axes * Two accelerometers are placed on each of the three 





Figure A-4 Configuration B 
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vehicular axes with the sensitive axis of one accelerometer along the 

axis in the positive sense and the sensitive axis of the other along 

the axis in the negative, -sense. The arrows in the figure indicate 

the sensitive axes of the accelerometers. 

Since A is at point a and oriented in the L direction, it 
V 11 

measures, from Eq. (A-9)s 


V"'‘V (0 v +0 ^ 

(A-14) 

By referring to Fig* A-3 and Eqs. (A«*7) through (A~12)* 
5' accelerometer readings can be determined * 

the remaining 

v ■ \ * * “V '* 

(A-15) 


(A-16) 

\ - «, * (4 * 4 > % 

(A-l?) 

\ \ ^ + 4 ) ! e 

V31 v 3 Vg v 1 

(A-18) 

\ “ \ + + 4 ) 

32 3 2 1 •- 

(A-19) 


These equations are now combined algebraically 



(A-20) 



i6. 


22 V 21 
+ *d 




(A-Sl) 


Then 


Thus 


■*7 " • ' 

x = . 
A e + A f 


“*1 + <,) 


< 


'*! + x 2 + x 3 




*1 " x ~ + x 


’g 3 i 

2 1 


*1 * *2 “ X 3 


*v, ■ < 


“*1 


+ x_ + x_ 1 / 2 


2 3 


) 


«V = 
2 


X 1 


l/2 

+ x 

_ — 2 ) 

*1 + X 2 “ X 3 1//£ 


(A-22) 

(A-23) 

(A-24) 

(A-25) 

(A- 26 ) 

CA-27) 

(A-28) 


These equations , .being wholly algebraic , are easy to evaluate 
on a digital computer. A difficulty arises, however, in evaluating B. 
since the square root of Eqs. J (A-23), '(A-2A) and (A-25) may have either ‘ 
a plus sign or a minus sign. This sign difficulty can be resolved through 
the use of auxiliary devices that may be less accurate and less costly 


than accelerometers . 



Then, from Eqs 


(A-22) 
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. (A-l4), (A-l 6 ), (A-l 8 )‘, (A-20), (A-2X), and 


\ = ^ + ^ "a ( A ‘ 2 9) 

\ “ ^ + X 2 *c (A-30) 


\ = \ + x 3 ^e 

v 3 31 J 


(A-31) 


If either the condition A a = S^, = £ a , and j? e = or 

£ a “ = £ c = J? a = A e = were satisfied, the equations presentet]. 

above would simplify. 

In configuration B, the accelerometers mounted parallel to the 
negative vehicular axes are Brought to the origin of the vehicular 
axis. as shown in Fig. A-k. In practice this is . impossible unless it 
were feasible to mount a three degree of freedom accelerometer at the 
origin, nevertheless, the configuration is analyzed here. The six 
accelerometers measure the following quantities: 

\l = \ " ( ^ 2 + ^3 } ** (A- 3 fe) 

Ay = Ay U-33) 

12 ! 

\ 21 = \ - <^ + *c 

Ay = Ay 

•22 v 2 


(A-35) 
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\ ■ \ + 4 > 

31 3 2 1 


32 v 3 


Since i and A,, are known [Eqs. (A-33), (A-35.) and 
1 v 2 3 
(A-3T ) 1 we can write 


2 2 

*v +s v 

2 v 3 


\ 

\ - *V 

,1 V ll- 


'd 2 

*V +S V 

V 1 v 3 


2 21 


2 2 
°T 

v 2 V 1 


~3 - V 31 


Then, as before, 

V = ( 


-x£ + x£ + x> 


1/2 


«V = ( 


+ x^ 


1/2 


<V = " < 

3 


x i + x 2 “ x 3 


1/2 


(a- 36) 

(A-3T) 

(A-38) 

(A-3?) 

(A-40) 

(A-41) 

(A-42)' 

(A-43) 


Configuration C utilizes 9 stationary accelerometers . It 
eliminates the ambiguity in the sign of the angular velocity S. Also 
no accelerometers are needed at the center of gravity. This method 
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gjjres the output as A ' A , a , o o and L' . Note that the 

’ / " 1 2 v 3' V 1 v 2 ' ‘ v 3 

mounting of the accelerometers must still 'satisfy the requirement 

that they he placed on the principal vehicular axis about the center 

of gravity. The mounting scheme is shown in Fig. A-5. 

Using Eqs . (A-7) through (A-12), we can write the outputs 

of these 9' accelerometers as (note that the alphabetic points nb 

longer coincide) 


\ 

12 

\ " (J V 3 ' \ *a 

(a-IA-1) 

\ 

32 

■ \ + \ + \ \ ) *b 

(a-44-2) 

^12 

= A Vi + 0^) *c 

(a-44-3) 

^32 

= a V s - - *v 2 «v 3 ) % 

(A-44-4) 

X 3 

= A ?i + ^ fV 3 ) *e 

(AJA-5) 

X 3 

X " ^ ‘ ^ *f 

(A-W-6) 

Xl 

^2 ‘ ( “Vg " ^ ^ 

(AJA-7) 

Xi 

■ \ + X 3 + \ V ^ 

(A-W-8) 

Xi 

= \g - x + \ %> ** 

(A-44-9) 
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Letting = % and aiding (A-kk-S) and (A-1A-!)-), 


\ + 

32 




32 _ 


V 


(A-45) 


Letting i = i and adding (A-it-4-l) and (A-W--3), 




2 


V 


(A-46) 


Letting ^ ^ and adding (A-Vi-8) and (A-44-9), 



V 


(A-4?) 


Three summations have yielded linear acceleration*. 

Letting and adding (A«44- 4) and (A-44- 6) yields % 

% + ^23 = \ X X % ■ X % + \ % " d 

= VV 2 V a 


* + '\' + ^V 

_ 32 23 v 3 2 

V 1 


2 ij, 


(A-k8) 
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Letting l = i and subtracting (A-1A-7) from {k-kh- 5) yields: 
6 g 




V \ *e + \ % *e + \ *e “ \ \ *e 

VV 2 V* 




Ay “Ay “ Ay + Ay 

v 13 31 1 3 

2 A„ 


(A-4-9) 


Letting and subtracting (A-Vt- 8) from (A-44- 1) yields:' 


\ 2 " \ " % + \ ^ X a “ \ 2 “ ^ \ g *a 

■ w 2 S* a 

■ \ -Ay -A_..-+\ 

fly = _1 IS IiS IS: CA-50) 

3 2 «a 


For Eqs. (A-45) through (A-50), the following equalities 


apply 


^c “ *1^ 

X d = f ‘f = 



A, 


(A-51) 



23 .- 


Then Fig. A-5 evolves as shown in Fig. A -6. Note that the 
arrangement is such that no .two accelerometers must overlap at a point 
(which is of course a physical impossibility) . 

Configuration D is of interest in that it utilizes eight 
accelerometers as shown in Fig. A-7- The equations governing these 
outputs are: 


31 

- \ 3 - \ - x \ 

(A-52-1) 

31 

= A^ * ^ 

IA-52-2) 

\ 2 

= \ + (Oy + *V *7 5 *e 

v 3 1 2 3 

(A-52-3) 

% 

= A V 3 ” ^ + *V 2 % 

(A-52-4) 

\s 

“ V + (®f + *V ) *e 

1 2 1 3 

(A-52-5) 

% 

= A v 1 “ ^*V 2 + f V 1 *V 3 ^ z f 

(A-52 - 6 ) 

\ 3 

= \ 2 - - <v 2 <V > * g 

(A-52-7) 


23 2 V 1 v 2 v 3 n 


(A-52-8) 




25V 



figure A-7 Configuration D 


Combining algebraically 
if Z e = t t = or; 


\ + Sr 



13 13 

2 


(A-53) 

if *g=\ “ 



+ 

v 23 23 

2 

- \ 

(A-5^) 

if * c = % = r 



\ +A v 
v 32 32 

- = \ 

(A-55) 


2 

also, let i - = A 


Then we can write the following equations using (A-52-l) 


through (A-52-8) 


= o - o 0 = x" (A-56) 

2& 2 V 1 V 3 



(A-57) 
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\ 3 I ^13 

2 (y 


% + \ \ ^ 


' V 

" " \ = X 4 


28 


x" + x," • V„ ' A V on ' V 


Then 


"7, 


•2 2jL 32 32 v 2 3 23 

2 ■ Hy 55 


xi' + x " \ \ - ■ K 

1 3 ■ *31 31 . T 13 T 13 

_____ - — T5— - 


Now from (A-56) 


% ‘ *1 

«V = _2 1 

3 


V 


• H - \ 

0, - - 31 31 

2 - 2A 

V 


Note that the angular information is 2^,2^ and 2^ . 

Equation (A-62) requires the integration of 2. [integration of 

V 1 

Equation (A-60)] before it can be solved. 
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(A-58) 


(A-59) 


(A-60) 


(A-61) 


(A-62) 
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APPENDIX B 

THE EFFECT OF LOCATING THE SENSORS 
AT A POINT OTHER THAN THE VEHICLE'S ’CENTER OF MASS. 

Alfred R. Schuler 

It is not usually convenient to locate the -sensors of^an 
inertial .navigation system 'at the vehicle's center of mass. . This 
appendix considers the effect, in an all-accelerometer system, 
of placing the ’ accelerometers symmetrically about’ the origin of-a 

. ■ . i 

coordinate system that is remote from and moving with respect to th 
center of mass of the vehicle. 

The equation governing the inertial acceleration at__an 
arbitrary point P accelerating with! respect to a 'moving’ reference 
frame has already been given by Eg.. (A-2). Restated for convenience, 
it is* 

Aj. = .Rj + fi x r y + ? v + 2fl x ? v + t x (ft x Jy) (B-l) ■ 

where 

A-j. = inertial acceleration of arbitrary point P 

"Ej = linear acceleration of origin -of vehicular coordinate. 

sysuem with respect to inertial space 
fl = angular velocity of vehicular system 
"r v = vector from vehicular origin to point P. 

Figure B-l shows the physical picture. The following two 
assumptions are made: 

1. The [m] axes, about which the accelerometers are plabed. 
remains in an orientation on- the vehicle that is paralle.- 
to the vehicle’s major axes. The vehicular axes [v]. 




30 . 


however j translates as a known function of time. Hote 
that [V 3 and [m3 always remain parallel to their respective 
coordinate axes . 

2. r..(t) = I r (t) + 1 r„ (t) + k r . (t) is known for all t. 

V 1 ' '2 V 3 

r^.(t) is stored in the navigation computer prior to 
launching or sensing devices are used to permit its 
calculation while in orbit „ 

To resolve Eq.. (B-l) into components we writes- 


R I = i V 1 +A v' 3 +A v- k 

V 1 2 3 

a = si i + si j + a. k 

V 1 2 3 

= ^ r 10 + + ( r 20 + ^ + ^ r 30 + 0 ^ 
*V = ^lO 1 + W + * 30 k 

K = ^lO 1 + + ^ 30 k 


where r-^, r 20’ r 30 are the of the ori Sin of the' [m] 

coordinate system with respect to the origin (center of mass) of the 
vehicular coordinate system, p, 1\ and £ are the magnitudes of 
the three coordinates defining a point in space with respect, to the 
[m] coordinate system and are assumed to be constants.. See Fig. B~2.. 
Equation (B-l) can now be expressed in the following fox-m: 
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Then 

1 2 v 3 V 1 v 2 V 3 

t(r 10 t |i)I + (r 2Q + T|)J + (r 3Q + £)k] 

+ ^10^ + ' r 20^ + 

+ 2(fly 1 + fi^T + ^k) X (r 10 i + ^ 20 J+ ^k) 

+ (<L I + J + SL k)' x. [(SL. 1 + SL J + fl k) x 

1 2 3 - X V 2 V 3 

[(r 10 + (i)I + (r 2Q + H)J +(r 3Q + £)k]} 

(B-2) 

Evaluating the three cross products in Eg.. (B-2) and collecting 
terms gives:' 

■*1 = + r 10 + % J- r 30 + ^ 20 + + 

' ' %/20 + X\( r 20 + H) + ^3^30 + 0 
- + ^ > ^ r 10 + 

+ j" [A v 2 + *» + ®v 3 (r io + *»> * \ (r 3 o' + C );+ \ 3 *10 

' ^ r 30 + W^IO + ^ + %\^30 + ^ 

-• (ft^+ nf^ ) (r 20 + 11)3 

+ k[A^ + r 30 + ^(rgQ + n) - 0^ (r 10 + p) 

+ Q V i r 20 " ^ ^10 + J V 1 J V 3 ^ r 10 + ^ 

+ W ^20 + (r 30 + e)] 

(B : 3)- 
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Equation (B-3) may Be mechanized By several methods, two of 
which are discussed here. First, an accelerometer arrangement identical' 
to 'configuration C described in Appendix A is utilized . Here, nine 
accelerometers are placed about the [m] axis. Mote from Fig. B-3 
that no two accelerometers are' located at the same physical point . 

By inspection of Fig. B-3 and Eq. (B-3), the readings of the nine 
accelerometers are: 

(Note that points and refer to the distances of the 

various accelerometers from the origin, of [m]) . 

to 21 \ 2 + r 20 + ( r 10 + ^^^3 + 

+ (fly^ ^ r 30 + J V 3 r 3U 

fgQ (B-4) 

A ’ m 21 k V 2 + r 20 + ( r 10 ■ + ^-^2 + ” ^ 1 ^ 3 ° 

+ ^10 - V^30 " + ^ 3 ) r 20 ^ 

*”31 = A V 3 + r 30 + ^^3 " ^ 10 + ^ 3 ^ + ^ 2 fl V ' 3 + ^^20 

(fly^ + fiy^) r 3Q + J2y^ ^ ^ 

*”l3 = \ + ^10 + <^3 + %> (r 30 + % 3 ) 

+ ft V 2 r 30 " n v 3 r 20 + T 20 

- <4 a + ^ 3 )r 10 


(B-7) 





3 1 ? 

*”23 = ' V , + r 20 + ^10 ‘ \ (r 3 0 f 5 + \*10 - \ f 30 - 

+ *^1^2 10 + \\ ^ 30 + ^2 ” ^1 + ! 1 

= V, + r 20 + + t V 1 < V 2 ^ r 10 + S-j " S^^O + 

+ - < V i P 30__" ( J V l + f V 2 ^ r 20 (B-8) 

*“12 + r 10 + S g r 30 " ®V 3 ' r 20 + ^ ^Vg^O ” fl Vj r 20 

+ \\ ( r 20 + + ^3 - (^ 2 + "lO 

\j_ + r 10 + ^ g + %^V^ T 30 + ^^2 * ^J^ 1 ’ 20 + ^ 

+ % ; 30 - % "20 - <<, + ) "10 ( B “ 9 > 

A ’“l2 = + r 10. + + t V 1 0 V'3^ r 30 + ~ ^ 3 ^ r 20 " ^ ^ 

+ V» ' • V 2 ° ' (fl * 0 + } r io ■ (B " 10) 

*“32 *V 3 + r 30 + 2 V 1 ^ r 20 + % 2 ) *V g r 10 + ~ ‘Vg^O 

+ n V i n V- 3 r 10 + ^2^3 ^20 + ^Kg * ‘ ^ + ^^30 

= *v 3 + r 30 + ^ + a v g n v 3 ^ r 20. + % g ^^^3 ‘ Sr 2 ^ r io 

‘V.^SO ' %'*10 ‘ ^ + ^ r 30 


(B-ll) 
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A ' m 32 = A T 3 + r 30 + ^ + ^20 )+ ^ £ V 1 S V 3 ‘ J V 3 ^ r 10 


+ n v , 1 r ao ‘ ^v/io “ ^V, + *VJ 


■30 


(B-12) 


X 2 

Then combining Egs.. (B-5) through (B-12), the linear acceleration 
and angular velocities are obtained. 


Affigi ” A m 2i 


8 \ 


«V + *V ‘V • “ ®1 

3 1 2 


(B-i3) 


A ff in. 


Ami 2 - A 


^-= = 


1 2 3 


(B-lli) 


^2 ,~ A ' m 32 jL + SL JL =9 (B-15) 

2% 1 2 v 3 3 

2 

*“21 “ *“23 = ^ ^3 + " *K g (fl V e fi V 3 " = ^1 

(B-i6) 

*”13'.’ *”12 = % 3 ’ ^3) ‘ \ \% - = h 

(B-17.) 

*”31 - *”32 " ^ ^3 " - %\ + W * = h 


How in Eg,. (B-l6) substitute Eg. (B-13) for 0£. 


(B-18) 
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“ \ e i ' \ (S V 2 V - 


~ K + \ 9 i 


% % " \ = e 4 


In Eg.. (B-17) substitute Eg. (B-l®) for 0_ 


*B.“ %3 ( \ V V ' 0£ \ 


02 + e 2 


<V *v + <V ■ = ®5 

1 3 2 5 


In Eg. (B-l8) substitute, Eg. (B-15.) for 9, 


%- 3 (S V^ % " \\ % ®3 


~^3 + % - 9 3 


W\ = 


The magnitudes of the components of Qy are given by 


6 1 - 0 2 
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e 3- ' S 4 




(B-23) 


e 5 * s 6 




(b-24) 


To obtain the components A^ , and , define the following 

1 V 2 V 3 

quantities ; 


2 2 

V + v 


= X 0 




+ c 


x i 


(B-25) 


< + 




= 


Note that these. fl's can he evaluated immediately after Eqs. (B-22)., 
(B-23) and (B-24) have been integrated. 

Combining Eqs . (B-4) and (B-5)" gives: 




Am 21 ' + A ' m 21 


r 20 " ®1 r 10 ‘ \ r 30 + X 3 r 20 


V T 10 + \ 


30 


IB -2 6) 


Similarly 


\ 


Am 12 + A ’“l2 


~ r l0 ■ e 5 r 30 ■ e 2 r 20 “ % s r 30 


+ y r 2o + x i r io 


(Br27) 
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and 


Am,„ + A'm-g 

\ 3 ~ “ 2 " r 30 ' ®3 r 20 ®6 r 10 


*V^ r 20 + % & T 10 + X 2 r 30 


If r, 1Q (t) = r 2Q (t) = r 3Q (t) 


= 0, then: 


\ = 


*“12 * A,m 12 


V 


Am 21 + A'sig^ 


V 


^32 + A ' m 32 


(B-28) 

(B-2 9) 
(B-30) 
(B-31) 


Thus, this reduces directly to the situation of sensors pla 
atiout the origin. 

It is noted that although the navigation computer must perform 
somewhat more arithmetic, there are only 3 integrations to he performed. 
There are no multiplications that. the computer must evaluate- (except for 
3 squares); Hote that if in general: ' 

r i0 (t) + 0 



r 30 (t) i 


0 
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but 


V 10 (t) = T 20 (t) = V 30 (t) = ° 


a 10 (t) = a 2Q (t) = a, ft (t) = 0 


30' 


Then 


V 


Am ie; +A ’ m i2 


®5 r 30 “ ®2 r 20 + X X r 10 


*"81 + A ' m 21 


"®1 r 10'“ r 30 + x 3 r 20 


Am, p + A''m_ e 

*V 3 = 2 _ ' “ S 3 r 20 ' ®6 r 10 + X 2 r 30 


(B-32) . 

(B-33) 

(b-34) 


This problem can also be approached through the utilization 
of six accelerometers that are placed with their sensitive axes 
parallel. to the [m] coordinate frame. This arrangement is shown 
in Fig. 4. -Wow using Eg. . (B-3) and Fig. B-t- we can write down by- 
inspection the .guantities measured with the six accelerometers: 

^11 A V 1 + r 10 + ®V 2 r 30 ' “v 3 r 20 + S 2 r 30 " %^ r 20 

W“ + W 30 ' ^2 + ^3 (riQ + " a) 


(B-35) 




*“12 = \ + '^10 + \ r 30 - \ r 20 + \* 3 0 - \ 


20 


W ® + W 30 ( ^ 2 + ^3 5 (r 10 - V 

Then subtracting Eq. (B-36') from (B-35 ) yieias: 


"'ll ‘ **12 - -<»? 2 * *. • <“? 2 * 4 3 > <1 


“*^1 - < .* - S 

2 v 3 


V*b 


Am, 


21 


■Vg + r 20 + “v^io " ^^30 + n v‘ 3 r io " 

+ W 10 + %V 30 " (s ^i + “V Crso + * c) 


*”22 \ 2 + r 20 + *V 3 r lO ' ^^30 + ^V^IO ' ^^30 

+ W 10 + % fi V 3 r 3° " + ^ r 20 ‘ 


Then Eq. (B-38) minus (B-39) gives: 


*“22 " *“21 “ + ®V ^ ^*4 + £ c^ 


to 22 “ *“21 2 2 

"W = vv x * 


in; 


(b-36) 


(B-37) 


(b- 38 ) 


(B-39) 


(b-4o) 



*®3i + r 3o + 20 ' *V 2 r io * 1 V 1 r ao " < V g r io 

+ \ J V 3 r lO + W*Z0 - ^ + <,> ^30 + V 

*“32 \ 3 + r 30 + Vj^O ' %*10 + S'^ao " ®V 2 r 10 

+ + ^ft/ 20 " + ^ 2 ^ ^ 3 ° " ^ 

Therefore Eq. (B-42) minus Eq. (B-4l) yields: 

to 32 ' *“31 = (n v i + U f + V 

or 


•*“32 ■ *“31 2 2 

“W^ = V V x i 


Theft hsing Eqs. (B-37)> (B-40) and (B-4-3) we have 


\ = 


X^ + X£ + X^\ 1 / 2 


' X 1 ■ X 2 + x 3 ' 1//£ 




+ x,; - x^ 


42.' 


(B-4l) 


(B-42) 


Ce-43) 


(B-44) 


(B-45) 


1/2 


(b-46) 



We can now write Aj in terms of known quantities. 

From Eq. (B-4l) 

\ 3 = *“31 ' r 30 " ^^20 + % 2 r X 0 + (% s ~ J V - L { V 3 ' r 10 

- \ + %^ r 20 + X 3 [r 3 0 + *e ] CB-47) 

Frpm Eq. (B- 38 ) 

AVg = *“21 " r 20 ‘ n V 3 r 10 + + W^IO 

+ ~ ^2^3^30 + X 2 ^ r 20 + A e^ (B-48) 

From Eq. (B-35) 

■\ 1 = A™!! - r 10 - Sg^O + " t®V 2 + Vv^SO 

+ - \\) r 20 + X i ^10 + * a 3 

Thus, the equations giving the .vehicles position; are easily 
mechanized when the accelerometers are placed about axes fm( narallel 
to the vehicle axes hut displaced from the center of gravity. Any 
change, due to fuel consumption, in the distance oetwesn. the sensors 
and the vehicle's center of gravity^ is. readily programmed -into the' 
navigation computer. 
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1A. 


ROTATING ACCELEROMETERS TO MEASURE LINEAR ACCELERATIONS 
AND ANGULAR VELOCITIES 

Alfred R. Schuler 


In this method for determining linear accelerations and 
angular velocities, each of two rotating rings carries an accelero- 
meter. The two rings are concentric and normal to each other, one 
rotating about the' V^ axis and the other about the V g axis as shown 
in Fig.(C-l)-. It is assumed that the rings are placed at the center 
of gravity and that each accelerometer is mounted with its sensitive 
axis perpendicular both to the axis of rotation and to a radial line 
from the axis of rotation to the accelerometer. 

Let the velocity of ring 1 be o)^, that of ring 2 be t»g . To 
evaluate the velocities and accelerations, first find the acceleration 
of a point D on ring 1 located at a rauxus y, rl irurn uue urxgxu ox uue 
vehicular system. The radius vector from .the origin to point D is 
given by 


' A d ~ ^dg *V + % 3 \ 

where A, and i, are functions of time and can be written 
d 2 d 3 

t d = A d COS (Djt 


(C-l) 


;c-2) 


A d sin aj-jt 



- 45 . 


v 3 



Figure .C-l ^Diagram .Showing -.the Position .Of the 3Srb Rotating 
; Bings :and -.the jLocation ;of the .’Acaeleromeiers .A- 
-:and Ag - L 



The 

coordinate 


velocity of the point D with respect to the vehicular 
system is given hy 


( 


d Ja 

dt 


) = 3, x £, 

v x a 


(c-3) 


and the acceleration with respect to the vehicular coordinate 
system is 

,2 7 

d 

( — p 2 - ) = (B, x (3) x jg, ) 

dt V x j. a 


(C-4) 


Expressing (C-3) and (C-4) in terms of their projections along 
the vehicular axis. 


^ dt ^ V®1 X ^ *d g + *V 

= "jy (»]_ ^ + Sy 

= -jy (U-j^ £ d sin ou^t + ky a> 1 cos m-jt (C-5) 

a 2 J, 

( S. ) _ _ _ 

dt 2 V = *V “l x ^“'V ®1 ^ + “l 

= - jy 0)! - \ !»! 

~ Q — v 

= -jy J8 d cos (D.jt - ty Sin (D-jt (C-6) 


It is shown hy Page that the acceleration measured in a 
moving system (the vehicular coordinate system [V 3 ) with respect to 


a fixed inertial frame is 
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(A^ = (Ejj + ( r) y + fl x (r) y + 2fl x (r) y + 12 x [fl x (r) y ] 

(0-7) 


Here, R is the distance from, the origin 'of the inertial frame 
to the origin of the vehicular system, r is the vector distance from 
the, origin of the vehicular system to the point at which the accelera- 
tion is being measured, and fl is the angular velocity of the vehicular 
system. Using Eg.. (C-7) to express the acceleration of point D in 
inertial coordinates. 


(A) 


I point D 



d 2 i. 


at 


+ 2a x ( 



) 

V 


) + a x U ). 

v a Y 

+ a x [a x (li ) ] 


(o-8) 


Expressing Eg. (C-8) as projections along the vehicular axis and using 
Egs. (C-5) and (C-6), we obtain the following expression for the 
inertial acceleration of point D: 


(A) I point D = V [ \ + \ *d Sin V " \ % C0S ^ 


+ 2(J) :L (fly cos iB^t .+ fly sin (Bjt) 


+ Ay (fly COS + fly Slltt 

12 3 


t)] 


+ jy [Ay - COS (Bj^t, - fly ^ Sltt (B ^ 



k&: 


an. 


i \ % eos «!* + p -v 2 *V % sin “i* 


(si|^ + ) i a cos Wjt] 


hy [Ay - sin o)jt + Sly eos ro-jt 

3 1 


a»i Sly S, & sin (Bjt + Sly Sly £ a eos o^t 


2 2 

(Sly + S^ ) SID BJjt] 


(o-9) 


Ay , and i represent the linear acceleration of the origin of the 
1 2 v 3 

[V] axis with respect to the inertial reference frame but resolved 
along the vehicular [V ] axes . 

Now consider the accelerometer Ag mounted on ring 2 which is rotating 
with angular velocity ujg. It is a radial distance J8 g from the origin 
of the vehicular axes . has projections on the and Vj axes and 
can be expressed as: 


As Ar A^ * Af ^feg 

i and % are functions of time and are given by: 
e l e 3 


^ = ^eosfflgt 


= jg g sin ffie t 


(C-10) 


(evil) 


The acceleration of a point E on this ring is: 

(A) I point E = h l \ - 4 As cos *2* + \ *e Sin u) 2 t 


SO 2 Sly COS 0)gt + Sly SZy^ £ & Sjjl ® g t 



" ^ n 7 s + -. E V^ ° os 

+ +(0^ < e cos m g t A e sin tog*) 

+ a« 2 (S2y^ jg e sin tDgt + JS^ je e cos togt) 

+ & m sia ®S fc + n v 1 eos “e 1 ^ 

+ *V -V, “ “4 -*e' sln u) 2 t “ V cos “a -6 

J 2 

" .V sin * 2 ^ + 'SV 1 fly < e cos »gt 
^ 3 

” C f 4 i + *§') & e sin ajgt] 


**■ 9 *' 


(C-12) 


Accelerometer A^ on ring 1 is assumed to have its* sensitive 

axis parallel to Vg at t = 0._ Shus,. it yill be directed along .-V p at 

t _ = > along " y 3 at t = ~~ , and along +Vg at t = =|H- , etc. 

In other words, the accelerometer lines up periodically along the 

posit! ve\and negative Vg and axes. Therefore, for t = 0. 

.Sff ^TT 2HTT _ 

«j_ *1 ? n = 1 5 2, ...... ® only the Sy component 

of Eg. . (C-9) is measured with accelerometer 
Also 

sia ( (U 1 ~~“™J * sin 2qtt = 0 
and 

cos |cn_ -i2lL_^ =, 


cos 2ntr » 1 
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Therefore, if we designate the quantity measured "by A., under 
these circumstances as we have: 

A n = \ 3 + \ + S 2 \ 3 A a ( c -«) 


For t 



5_ rr 9 ft 
2 3 2 



n = 


only the component of Eq. (C-9) exists and this is in a negative 
sense. 


Also 


and 


"‘■fc 1 

~ (»i » 


Designating- the quantity measured hy A 1 for the cohdition 
that t = ^ n „~ 1 ' to he A,_ we have: 


2lu-| 


12 


A 12 = ‘ %' + ^v” %■* 


_a v 2 + 9 v 1 % ' % % % 


(C-l!f) 


For t = 


- _IL_ 3rr 5w 


(2n-t 1 )tt 


h 2, 3 


again only the component of Eq. (C-9) exists and it too is in a 
negative sense. Call this quantity A-^* 
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f ^)- » 

4i 


Therefore 


- [ \ 3 *a - \\ % ] = ”\ 3 + ^ * a + % % % 


F° r t = -|£- , , ...... ^ i £g ~ £ ^ , only the J y component 

. of Eq . (c-9) is -detected lay accelerometer Designating this quantity 


as A n » and noting that 


“-K- lis ^r k -) ■ - l 


we then have ; 


^1 = + V % ‘ V 


We now follow the same procedure for accelerometer at 
point E on ring 2 as we- did for accelerometer A^ at point D on ring 1. 

Thus, for t = . -we note that only the £* r component is 

' - -V 1 

detected with Ag which from vu -12) is: 


*21 = % - %«e + \ % *e 


(c-17) 
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Also, for t 


(4a +- 1 )tt 


only the minus l^. qomponent 


detected . 


Designating it ty A^', Eg.. (C-12) yields,: 


(C-18) 

From 

( 0 - 19 ) 


(C-20) 

Assume that and also that the ring frequencies' and 

(Mg are equal to u) so that the sampling times are identical. Bren 
repeating the eight equations through Ag^ for ease of manipulation: 


A 11 

= \ 
v 3 


% 

+ % 

% 

% 

(c-21) 

A 12 



% 

- v 
v 2 

% 

%’ 

(c-22) 

A 13 

= ' V 3 

+ \ 

% 

+ 

V 2 

S 

% 

(c-23) 

a i 4 

= s 


A d 

"S 

% 

A a 

(C- 24 ) 

A 21 

= % 

-% 


+ t v. 


5 % 

(C- 25 ) 


A 22 " ~ A V 1 " “Vg ‘ % 1 % 3 l e 
only the minus component is read. 

A 23 = "V, " Sg A e + ^ (V, l e 


rf- 3 }rr 

For t = ~ 2 ^ — , tne component Ag^ is detected hy Ag: 


Fort = i^Lik 
®2 

Eq. (C-12): 


& 2k 


\ " \ A e - \ \ 
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A 22 = A V 1 " %• ! V 1 % 3 % 

A 23 = _A v 3 " % + S V 1 % 

k 2h = \ 1 * K 0 1 <l ” ^ ®v 3 A a 

Subtracting (C-23) from (C-21) and dividing by two yields: 



Likewise, working with other of the above equations : 


A 24 ” A 22 


“ V 


A lk " A 12 


V 


A 11 + A 12 + a 13 + A ih 

-TT— - V 


~ A 21 - A 22 " A 23 ' A 2k _ • 

ns = v 


To obtain SI , (c-21) and (c-23) are added 
> V 3 


V + A 13 


\ + % % 


(C-26) 

(C-27) 

(C-28) 

(C-29) 

(c-30) 

(C- 31 ) 

(C-3S) 

(c-33) 



(c-34) 


Then solving for 


V 


A 11 * A 13 
2 Jfc, 




V 


It is seen that JL' i A, < i- » fl. ah'! 4, are found directly. 

.1 2 . v 3 1 2 

,A single integration of fL. and sL yield the angular rates SI and Sl^ . 

1 V 2 V 1 V 2 
SU , given by Eq. (0-3^) , must be found by solving a simple algebraic 
3 

equation . 

This method has as its primary advantage over fixed accelerometers 
the elimination of at least four additional accelerometers and associated 
electronic gear (amplifiers and torquers). It -would be useful if the 
following conditions are satisfied: 

a) Linear acceleration and angular velocity are slowly 
varying functions of time so that to is sufficiently 
-.low to insure that the time lag in the accelerometer 
output (due to the restriction on frequency response 
of the accelerometer) does not introduce errors. 

•b) Sufficient power is available to drive the rings on which 
the accelerometers are mounted. 

c) Practical sampling circuits are not too difficult to 
implement . 

d) Packaging presents no major difficulties. 

REFERENCE 

1. Page, Leigh: "Introduction to Theoretical Physics," 

D. Van Nostrand and Company, 19Vf. 



55 . 


APPENDIX D 

INPUT-OUTPUT ERROR RELATIONS 
OF Hffi GRAVITY COMPUTER 

Anthony J. Grammatieos 

The function of the gravity computer Is to generate the 
gravity acceleration based on the vehicle position information 
available from the navigation system and the universal law of 
gravitation which the computer simulates . 

The position, however, is not known precisely and therefore 
position errors enter the computer. As a result, the computed gravity 
components contain errors . The output gravity errors of the computer 
are related to the input position errors. This relationship is needed 
in studying the stability of the navigation system. The present 
appendix establishes this relationship. 

Let nu, x^, y^, z^ (i = 1, 2, ... n) be the masses and 
coordinates of the centers of mass Of the celestial bodies contributing 
to the gravity acceleration at the true position x, y, z of the 
vehicle. Then the true components of the gravity at this point are: 

(D-la) 


V - 
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g = „ S ym . . — - 
z i=l )’ R| 


(D'lc) 


where 


R x = [ (x-Xj) 2 + (y-yj) 2 + (b-Zj^) 2 ) 


l/2 


(D-2) 


and 

_8 0J31^ 

•y is the gravitational constant , 6.67 x 10 — — g 

grosec 

Since the true position x, y, z of the vehicle is not available, 
the gravitation acceleration is being calculated on the basis of the 
approximate position information Sy , s^_, Then the calculated 

gravity components are; 


T, y m 
i=l 1 


s'? 

X 


(D-3a) 


g = - 2, v nu 

V iix ^ i 


I S y - y i 


S 3 

1 


(D-3t) 


- - S *y m„ 

i=l 1 


(D-3e) 


where 


s i = ^ ^ s yf x i) d + ^ s y" y i^ + ^ s a" z i^ ■* 


^i 1 / 2 


(D-^) 


Let the position and the gravitation errors hes 


Sx = s x ~x 


8 y = s y -y 


8 z = s z -z (D-5) 
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6s x = ^ C -S x 


5s y 


6g, = g -g„ 


(35-6) 


where the g's are given hy (D-l) and (D-3). 

How expand Eq. (D-3) in Taylor series about the point (x, y, z) . 
Take for example Eq. (D-3a). 


g (s , s , s ) 
xe x y z' 


s y = y 




(s x - x) 


s - y 
y J 


3 

3s 


(s r - y) + 


Sy = y 


3d 


(s - z) 
v z ' 


s y = y 


+ higher order terms in (s -x), (s -y), (s ~z) 
x . y z 
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* s x Ss x 


- S- y m. 
i=l 1 


. 1=1 V 1 BSx ' + iwf + (■ ^f *' 2 


= r L. y ffii 


t(s x - Xi ) + (B -y ± ) 2 + (s z - z.ff 


( V X l) I t(B x ^) g * (s y -y.) 2 V (Sg-B^ 2 ] 2( V x.) 

* ^y- y i )2 + ( s z - z i ) 2] ' 6/2 


n 1 . 3(s -x ) c 

s vm . s? 3""^ 

i=l l s. 

l 


= - S y m. 


3(V x i)‘ 


3(x - x ± ) c 
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Similarly 


3-§xc 



n 

E 

1=1 


V m i 



3(x-x ± ) (y-j ± ) 


(D-10) 




3s z 



- E 
i=l 



r - 3(x-x ± ) (z-z..) 



(D-ll) 


In view of Egs. (D-5)> (D-9)? (D-10) and (D-ll) and dropping higher 
order terms Eg.. (D-7) takes the form 


n v m. j 

y -A ’ 

1 (x - X. ) + 

1 1 

1 1 - 

3(x - x..) I 2 - 

i. O i 

i=l Rf 


R? 

x | 

l. 

_ 

i J 

3(x - x ) (y - y. ) 

3(x - 

x ± ) (z - Z ± ) 

R i 





6x 


6z> 


I (D-12) 


Ihen from Eg. (D-la), (D-12) and (D-6) we obtain 


n *y m. 

6Sx = 8xc ' Sx = ' hi .if 


i - 
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and by analogy/ 
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Then Eq. (D-13) can be rewritten as follows; 
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5s x = - ti x 8x + v 8y + p Sz 
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(D-l8b) 


In most cases of an interplanetary flight the vehicle will 
"be 1 in the sphere of influence of the sun only and equations (D“13) 
to (D-18) will he simplified considerably * Whenever* several 
celestial bodies significantly contribute to the total gravitational 
field in the vicinity of the vehicle* it is convenient to introduce the 
concept of an equivalent celestial body with an equivalent g/R ratio of 
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S 

i=l 


•V “i 



(D-19) 
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Note that an- error in the measured value of any position 
coordinate will cause an error in each of the components of gravity 
errors and that the coefficients of Eq. (D-l8) are timervarying . 

For'a unique center .of attraction (the two tody problem) 
the coefficients given "by Eqs. (D-14J, (B-15), (D-l 6) and (D-17) have 
.hounded values as follows: 


if"" s » V 

VI 

(D-20) 

3g 

2R V j p 3 T €. 

3g 
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In the case of’ a unique center of attraction some of the 
previous equations can he rewritten as follows: 


where 



‘ ^ 

(D-22a) 

“ ^ 

(D-22b) 
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,(D-22c) 
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using Eq. (D-23) in Eq. (D-l8) yields 



(D-2!t) 

A careful examination of the Taylor expansion employed previously 
reveals that the first term in Eq. (D-2^) represents the gravity errors 
due to the position errors included in the numerator of Eq. (D-3) 
and that the second part of Eq. (D=2k) represents the gravity errors 
due to the position errors included in the denominator of Eq. (D-3). 
Returning to Eq. (D-l8) we can write 

6g = G §R (D-25) 

where 



Equation (D-25) is the basic result of this appendix. 
Recapitulating, Eqs . (D-l) and (D-3) can he written as 

g = G x (x, y, z, x ± , y.., z ± , n^) (D-l) 




or in view of Eqs. (D-5) and (D-6) 


'6k . 


81 = g e -g = G 1 (x J y, z, 8x 3 8y, 8z, x ± , z ± , m ± ) 

“ y, z, x^ y ± , z ±3 n^) (D-27) 

Here ( . ) represents the functional form of the universal 
law of gravitation. 

Among others, Eq. (D-27) gives the dependence of 5g upon 
6x, 8y, 8z which is needed in the error analysis of the overall 
navigation system. Such a dependence, however, is nonlinear and cannot 
he used for an error analysis directly. Then what we dd is to find 
an approximation of Eq. (D-27) by expanding Eq. (D-27) in Taylor 
series in terms of 8x,' 8y, 6z about the point x, y, z, and retaining 
the linear terms only; this can be done based on the assumption that 
gx, 6y, 8z are small. 

This procedure leads, to Eq. (D-25) which is valid for the 
errors only . Apparently G ^ . The elements of G .Vary slowly with 

time and act as gains upon the elements of SB to generate the elements 
of 6g. 



MTEtlDIX E 
A- Grammaticos 

E.l THE MECJHAMIEATIOIf OF THE NAVIGATION EQUATIONS IK THE 
INERTIAL FRAME 

The problem of the mechanization of the navigation equations in^ 

different coordinate frames has been studied by Bodner and Seleznev" 1 " 

2 

and Krishnan in their works on gimballess navigation systems „ Of 
all the mechanizations that were studied, the mechanization in the 
inertial frame was found advantageous and recommended for further 
study. 

Figure E-l shews such a mechanization whose fundamentals 
are presented here. Assume that the space vehicle on an interplanetary 
mission is moving along a path around the sun. The vehicle is subject 
Jo' the gravitational forces and the thrust of its. engines. 

The origin of an inertial frame is placed at the center of 
the sun and its x, y, z axes are fixed with .respect to the fixed 
stars , A vehicular frame is rigidly attached to the vehicle and 
follows its motion. 

The relative orientation of the two frames can be established 
from on-board acceleration measurements made with body-mounted 
accelerometers . The inertial acceleration a measured along the vehicular 
axes is resolved along the inertial axes by means of a direction- 
cosine computer, yielding the components a , a , a. of the inertial 


acceleration . 




Figure System Mechanization 
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Let i(g x> g^j g z ) be the gravity acceleration at the position 
occupied by the vehicle. Then the vehicle is subject to a total 
acceleration 

a + g = R (E-l) 

where R(x, y, z) is the position vector referred to the inertial 
frame. Equation (E-l) expresses the Hewton law of inertia. 

Integrating both sides of Ea. (E-l) once yields the velocity 
R(x, y, z) and integrating twice yields the position R(x, y, z) of 
the vehicle provided that a and g are known, a can be obtained in the 
way described above j g is obtained from the universal law of gravitation 

g = ^ E (E-2) 

R 3 

which is simulated by the gravity computer; y is the gravitation constant 
and m is the attracting mass. By applying Eq. (E=2) we introduce the 
feedback loops shown in Fig. E-l. 

Finally, assumed initial conditions on velocity and position 
are shown in Fig. E-l. 

Recapitulating the accelerometers measure components along the 
vehicular axes of the inertial acceleration due to the engine thrust. 

These components are transformed into components a x , a^, a z along 
the inertial frame. 

The sum of the thrust acceleration components a.^, a y , a z and 
the corresponding gravity components g x , g^, g z must he equal to the 
acceleration components x, y, z of the vehicle. 



Hie vehicle acceleration components x, y, z are integrated once 
to give vehicle velocity components x, y, zj proper initial conditions 
are introduced. 

She vehicle velocity components x, y, z are integrated again to 
give the vehicle position coordinates x, y, z; proper initial position 
coordinates are introduced again. 

Finally une current vehicle position coordinates are fed into 
the gravity computer which generates the gravity acceleration needed at 
the input of the first integrator. 

From a system point of view, this mechanization contains two 
distinct parts cascaded together. The direction-cosine computer is the 
first part; the other is what one usually calls the navigation loop. 

E .2 ERRORS "iN THE KWI6ATI0H SYSTEM 

Figure E-l shows the ideal navigation system. If the aeeelero-' 
meters measure acceleration without errors_and the initial alignment 
of the system is perfect then the system output is identical to the 
true position of the vehicle. 

In practice, however, it is impossible to make error-free 

acceleration measurements and perfect initial alignment; as a result 

errors are introduced into the system. Because of the errors 

6^ T (6 a v j_; & a v 2 J ® a v3^ in acceleration measurements the output of 

the direction-cosine computer contains- errors 6&(8a , 6a , 6a ) which 

x y z 

enter the navigation loopj . these acceleration errors together with the 
initial condition errors give rise .to position errors -6R(Sx, 6y, 6z); 
the position errors i enter the gravity computer giving rise to gravity 
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errors 81 ( 5^3 Sg^j 5g 2 ) which ip turn are fed taels into the system. 
Therefore, a staBility analysis of the system is necessary. 

Because of the cascaded connection Between the direction "-cosine 
computer and the navigation loop the stability analysis of the system 
can Be split into staBility analysis of the direction-cosine computer 
and staBility analysis of the navigation loop. 

Here we are concerned with the stability analysis of the 
navigation loop. 

Figure E-2 shews the navigation system and the corresponding 
signals; each signal is represented as the sum of its true value plus 
an error. By inspection of Fig. 2-2 we can write 

a + 8a + I + §g <=> R + 6fi (E-3) 

Using Eq. (E-l) in Eg.. (E-3) gives 
8a + 8§ - 6R 

or 65 - 6g = Saj ' ( E " 5 ^ 

Equation (e- 5) is the error equation for the navigation loop. 
MOtc that 6a, acts as a forcing function whereas 8l depends on SR, the 
dependence of 8g and SR was estaBlished by Ea.’s (0-18) or (TS-25) and 
(D-26) . We recall that; 


p- — 


— 



8S X 


V p 


6x 

8S y 
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V T 


6y 

^z 


P T ‘^Z 


$z 
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Assraning that p , p > p , vj p 5 T vary slowly with time, we 
x ■ y z 

ean Laplace transform Eq.. (E-7) 



(E-8) 

Now the stability of the navigation loop is determined by roots 


of the characteristic equation „ Consider first the determinant 
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2 

S -Hi 
r x 

“V 

"P 

D(s) 

= 

“V 

2 

s +11 

“T 



-p 

“T 

2 

S 4 j| ( 


(E-9) 


Expanding gives i 

D(s) = (s -tp x ) (s S +ii y ) (s 2 +^ z ) - vTp - vTp - p 2 (s 2 +p y ) -r 2 (s 2 +ji x ) -v 2 (s 2 +p 
or D(s) = s 6 + (p x +p y +p z ) + (p x By. + p y n z + |i z M x - v 2 - p 2 -.T 2 )s 2 

" (M x t 2 + p y p' 2 + \i z v 2 - p x |i y + 2 vpT ) (E-10) 

In view of Eq.'s (D-l4), (D-15), (D-l6), (D-17) and (D-19), 
we obtain: 

+ + '= 0 (E-ll) 

P x f* y + M y P a + H z |x x “ v 2 - p 2 -/ = -3 Sg- ( E - 12 ) 

R 

M x T 2 + My p 2 + n z V 2 - M x Hy n z + 2 vpT = 2 (E-13) 

R 

Also let | = m 2 (E-14) 

Then Eq. (E-10) can be rewritten as: 

D(s) = s 6 - 3 dig s 2 - 2 u>g (B-15) 

D(s) = (s 2 + ra 2 f (s 2 - 2 u) 2 ) 


or 


(Erl 6) 
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Let 
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s +u 
“x 

-v -p 
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-v 

s2+fi y “ T 


(E-17) 




L -p 

s Sf 



s'-inee |q1 

= B(s) ^ 0, the 

solution to Eq. (E-8) is 
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(E-18) 
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(E-19) 
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(E-SO) 

Equation (E-20) suggests that (the characteristic equation of the 
system (E-8) is 

(s 2 + a, 2 ) (s? - £b 2 ) = 0 (E-21) 

The roots of the characteristic equation are ±o‘u) g and ±V2~ra s . 
This means that the transient response of the navigation loop consists 
of a sinusoidal term of hounded amplitude due to the ±ju) poles and a 
hyperbolic cosine term, which increases with_time, due to the i/? ® s 
poles; therefore the system is unstable. Small errors in the initial 
conditions produce time-increasing error in the output, and the 
accuracy of the system deteriorates with time. 

In order to improve the accuracy of the system we must stabilize 
the system by damping the errors . Methods • of damping are presented in 
the following section. 
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E.3 damping the system errors by means of reference trajectory 
INFORMATION 


The reference trajectory is defined as the ideal • predetermined 
trajectory for a specific mission. During the execution of the. mission 
one does his test to keep the actual trajectory as close to the 
reference trajectory as possible. Hence, at a given time t the difference 
between the corresponding positions on the reference and the actual 
trajectories is small. 

The method of damping proposed here is based on the closeness 
of the two trajectories mentioned above and on the assumption that 
reference trajectory information is available on board the vehicle j 
this is usually the case since this information is needed for a number 
of purposes including guidance. Let 


R 


r 



be the position that the vehicle should 
occupy on the reference trajectory at 
time t. 


8R = R - R = 
r r 


8x„ 




6z_ 


and, anticipating its use. 


be the position difference between the 
reference and actual trajectories at 
time t. 

form the difference 


1 

~ 8x r 

- 8x 

R r - r = R + 5R r - (R + 6R) = 6IT. - 8R = 

6y r 

“ 


“ 8z r 

- 5Z- 


(E-22) 
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Figure E-3 shows the system mechanization of Fig. E-2 with 
reference trajectory damping added. In this mechanization, the 
reference trajectory position R r is compared with the calculated 
position s’ and their difference "5^ - s’ is modified by the compensation 
matrix 1. The difference R r - s' is given hy Eq. (E-22); the compensation 
matrix H is to be determined so as to achieve damping. 

How we can write the error equations by inspection, of Fig. E-3; 
they are 

R + 8R = I + 5l+a + 6a+S (SRj. - «R) (E-23) 

or „ 

SR = 6g + 8a + H UR r “ 8R^ (E-2^) 

since R = g + a 

Ey.^>-j-on (E-2 1 !-) can be rewritten 

8R + H SR - G 8R = 8a + H 8R r ■ (E-25) 

Equation (E-25) is the basic error equation of the system in 
Fig. E-3. The terms 8a + H 8R r constitute the forcing functions and' 
hence the error equation of the autonomous system is: 

6S + 3‘ SR.- G 8R = 0 (E-26) 

Equation- (E -26) determines the performance or the system. 

Careful inspection - of this equation aids- in determining a suitable 

form for 3 . The damping of the; system is directly related to the presence 

of- derivative terms tk 6R; so H must contain terms of the form k - 

dt 

or kp -if p = -ir . The speed of response is related to the coefficient of 
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8R (i.e. the undamped natural frequency) which we would like to be 
able to control as much as possible . Since we have no control over 
G we introduce a constant part K in H so as to obtain the term (K-G)8R 
whose coefficient K«G is adjustable . 

With this in mind let us suggest the following form for H: 


(2k x P +Kp 


yx 


L- k 


xy 


(2k j +k>) 


■k 


zy 


yz 


( 2 k v + k’)- 


(E-27) 


where all the k's are constants or slowly varying with time. 

Recall that G is given by Eg,. (D-26)j introduce both G* and IT 
into Eg. (E-s 6) take the Laplace transform and expand to obtain: 

(s 2 t2k x s + k^u x )8x + (k xy ^)8yt(k xz -p) 8Z =Xt8a x +(5k Sc p+ki)8x r+ k xy 8y r +k xz 8z r 3 

(E-28a) 

(k^-vJSx+Cs^yS+k^JSy+Ck^-rJSz =X E8a y +k yx 63£ r +(2k y p+kpSy r +k yz 8z :t 

^E-28b) 

(k zx -p)8x+(k zy -T)8y+(s 2 +2k z s+k'+ t i 2 )8z =/[8a y +k zx 6x r +k zy 8y r +(2k z p+k’)8z r ] 


(E-28c) 


The characteristic equation of the system (E-28),is: 


D(s) = 


s +2k s+k®+u 
x x p x 


k » v 
yx 


k - v 
xy 


s +2k s+k*+u 
y y *y 


k - t 

zy 


is. _ k 

XZ P 


k ~ t 

yz 


s +2k s+k® +ti 
z z 


= 0 


(E-29) 
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Since \ij p 5 t are slowly varying terms, it is possible to generate 
k's such that: 


k 

xy 

= k 
yx 

s v 

(E-30a) 

k 

yz 

u 

= T 

(E-30b) 

k 

zx 

~ k 
xz 

= P 

(E-30c) 


Then the characteristic equation becomes: 

(s 2 +2k s+k'+ji ) (s S +2k s+k'+jj ) (s 2 + 2 k s+k’+u ) = 0 (E-31) 

j J J z z z 

We' like the roots of the characteristic equation to be complex 
conjugate with negative real parts ; hence the following conditions must 
be satisfied. 


k x >° 

(E-32a) 

k x + ^x > 0 

(E-32b) 

1 

£ 

+ 

£ 

A 

o 

(E-32c) 

k > 0 

y 

(E-33a) 

k* + u >0 

y p y 

(E-336) 

k 2 - (k ! + (i ) < 0 

y y V 

(E-33c) 

k z >° 

(E-34a) 

k- + Mz > ° 

(E-3®) 

k z ‘ (k z + »* z ) < o 

(E-3^c) 



80.. 


To show that /these eonditions can he satisfied and lead to a 
desirable situation we recall the condition (D-20); that is: 

- f 6 * IV V * I {D ' 20) 

Let us look at the inequalities (E-32), keeping in mind that 
§ = u) 2 is the undamped frequency of the basic navigation loop without 

K S 

any external information, and that this quantity is very small for an 
interplanetary flight. 

Combining eonditions (E-32b) and (E-32c) we obtain 

0<k^'<k^ + sJ x (E-35) 

or 0 < k^ < k£ - (E-36) 

2 & 

given that min - - -§■ 

X K 

Clearly the condition (E-36) can be satisfied by a proper choice 

0f k x* 

Finally the conditions (E-33) and (E-3k) can be satisfied in a 
similar way. 

Then the transfer matrix H takes the form 



(E-37) 


where the constant gains satisfy the conditions set above. 
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Prom Equations (D-15)? (D-l6) and (D-17) note that the transfer 
gains v> pj T in Eq. (E-37) depend on the true coordinates of the 
vehicle x, y, z which are not available; instead we generate these 
gains from the output $ z of the system. As a result the 

generated gains are v+8vj p 4 8pj t + 8t and when these values are introduced 
into Eq. (E-38) the terms involving 6v s 6p, 6t are of higher order and 
can he dropped from the equations. 


How examine the possibility of avoiding the use of variable 
gains; a simple choice for H with this property might be the following: 


H = 


2k p+k* 
x^ x 


2k p+k' 

y * y 


j 


(E-38) 


Expanding Eq. (E-26) and taking Laplace transfbrms we obtain: 


Cs 2 + 2k x s + k^ x 3 8 x -v8y - p8z = Jf [sa x +(2k x P+k^)6x r ] (E-39a) 

-v8x + [s 2 +2kyS+ky+p y ]8y - t8z = JP [6a y +(2k y p+kp6y r ] (E-39b) 

-p8x - T 8y + [s 2 +2k z s+khi z 38z =£ [8a z +(2k z p+k^)8z r 3 (E-39e) 

The Characteristic equation of the system (E-39) is: 

[[s 2 + 2k x s + k^ x ] -v "Pi 


D(s) = 


[s +2k s+k’+u 3 

y y y 


[s +2k s+k'+tt 3 
z z z 


= 0 


(E-to) 
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Now let 

k x = k y = k, = en and = *£ = k* = J 2 2 (E4l) 

and expand Eg.. (E-it-O) as follows: 

(s 2 + 2 gns + fl 2 ) 3 + (p x + P y + v z ) (s S + SJJls + fi 2 ) 2 

+ H y + ^y |l Z •< Jig V X • V 2 ~ P £ “ T 2 ) (B 2 + 2C0S + f( 2 ; 

+ ^x >*y ^z ■ ^ “ **y f 2 " ^2 y2 ■ = 0 (E-te) 

Recall that the elements of G satisfy the following relations: 

+ Ry + W z ~ 0 (E-43a) 

'2 

M x M-y + J* y fi z + u z p x - V 2 - p 2 - T 2 = - ^§~ (E-43V) 

K o 

2 2 2 2g J / i \ 

U x *V ^ T “ P “ V ~ Sv P T = ~ ^ 3 “ (E-43c) 

where 

R = (x 2 + y 2 + z 2 ) 1 - 72 (E43d) 

g « (g x + 8 ^ + g z ) l/2 (E-43e) 

Hence Eg. (E-t2) can be rewritten as follows: 

2 

(s 2 + 2£fls + G 2 + (D 2 ) (s 2 + 2gGs + n 2 -2cb 2 ) = 0 (E-44) 

where a > 2 = g/R 

It is desirable that all the roots of Eg. (E-V)-) have negative 
real parts for stability reasons and are complex conjugate in order to 
achieve fast response. It is possible to satisfy both reguirements 


as follows: 
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Let s. , s 2 be the roots of 

2 2 2 
s + 2Jfls + fl +' o)g = 0 


and s^j s^ he the roots of 


2 2 2 

s + 2t;ns + a - 2 Wg = o 


s l 5 s 2 = _Cfl 4 ^ fl2 - n2 - »s ^ 


s 3> s i). = -£ n 4 t £ 2 n 2 - a 2 + ai>g]_ 


In order to satisfy the requirements on s^, s^, s^ and s^ 


we demand that 


Gfl > 0 


2 2 2 2 
C a - U - < 0 


2 2 2 2 

r a - a + a» < o 


How Eqs. (E-!+5b) and (E-4-5e) yield 


c -k -§ 

a 


2 “s 

C - 1 < - 2 

a £ 


or combining them 


2 2 

<» s % P 

< 2 -§ < 1 - C 

a a 


(E-45d) 



and 


(E-46a) 

(E-46b) 


All conditions (E-45) are satisfied if we set 

C > ° 

ti > 0 


1 - 


C 2 > 


0) • 

^ -f 

o 


(e-1i-6c) 


A typical value for the damping would he £ =0.7; then from 
(E-46e) we obtain 


1 - 0.7 


2 !k 

n 2 


or 

a > a j) s 

Note that the period T = — is of the order of several months 
s a) B 

in- the case of interplanetary flight; hence io g is small. A large fi 
will increase the undamped natural frequencies of the system thus 
decreasing the period of oscillations considerably . Such an effect- 
is desirable. 

From the above discussion it is clear that both the damping £ 
and the undamped frequencies can be adjusted with considerable freedom 
so that a desirable transient response can be achieved. 

Looking at the right hand side of Eq. (E-28) we notice that 
the extra forcing term H SB^ was introduced; this results in increased 
steady-state errors. Given, however, that is relatively small 

the increase in the steady-state error is small. 

This method of damping can be used in a time interval during 
which it is ■ impossible to receive external information for some reason 
(as, for example, in the case where an observed-celestial body is pbscured 
by another). The method hasjwo important features: it damps 
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the navigation system and is completely self-contained. Its accuracy 
depends on the closeness of the actual and reference trajectories 1 
during the time of operation in this mode. 

Finally, let us redraw Fig. E-3 for the sake of simplification; 
this is shown as Fig. E-4. Here, instead of injecting the error signals 
2£0s [ 8R r - 8R ] into the acceleration node of the system, we 'inject 
the error signals 2£fl [ SR^ - SR ] into the velocity node. Obviously 
this does not change the situation at least as far as the characteristic 
equation is concerned. 

The values of the gains shown in Fig. E-4 are: 

k = *; = k = « 2 

2k x = 2k y = 2k z = 2 C 0 



R*f S Rf 


r+Sr 
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e.4 damping the system errors by means op external velocity information 


The method of Doppler velocity information for aiding airborne 
inertial navigation systems is a relatively old one. It has given 
satisfactory results in the case of terrestial navigation. 

In the -case of inertial space navigation, the requirements are 
more severe and radar Doppler is inadequate . 

Efforts to cope with the difficulties involved have led to new ideas 
of which the most important in the area is the optical doppler method. 
Feasibility and accuracy studies of optical Doppler^ led to 
encouraging results; the test of a breadboard model verified these 
conclusions . 

It is here assumed that such a practical device will be available 
for use on board the veh’~le. It will measure components of the vehicle 
velocity in the direction of a fixed star. 

Now, define: 


dx 


ay 


dz 


the velocity of the vehicle as 
measured by Doppler devices. 




6y 


dx 


6V 


ay 


6v, 


dz 


the error in the Doppler velocity 
measurement 
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then 


T d = R + fiv d 


x + 8v, 


ax 


y + Sr 


dy 


+ 6V; 


dz 


We also have 


E + SR 


and 


v a ‘ 


x + 6x 
y + Sy 
z + Sz 


5v dx - 8x 
5v dy - sy 


8v dz - 6z 


since s = H + 8H 


(E-W) 


Ihe suggested system configuration is shram in Pig. E-J. 

She Doppler velocity measurement v^ is compared with tfie velocity I" 
obtained from the system. The errors 6v a - SR resulting from this 
comparison are modified by the compensation matrix H (as yet specified) 
and -Orem are fed into the acceleration node of the system. 

The elements of the matrix 5 will be selected so as to eliminate or 
reduce -the time increasing and oscillating errors in t|ie system. 

From Fig. E-5, the error e^mtioij of the system is: 

- SR = S SR ' + ga + H (5v d - SR) 


,(E-48) 



FROM 

ACCELEROMETERS 



S x *X+§X 

S y =y+8y 


S z =z+Sz 


Figure- E-5 Velocity Damping 
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or 


8R,+ H 8R - G 6R = 8a + E 8v, 


[E-4-9) 


'compare Eq. (E -49). with' Eq. (E-25) of the previous section and note 
the- similarity -between them.. Let us choose a similar form for H. 
'That is, let 


(2k x p) 


/yx 


k 

■ xy 

k’ 

(-2k r + -f ) 


zy 


yz 

k’ 

+ ;r ) 


(E-50) 


' where all the k 1 s are constants or slowly varying .with time . 

Introducing, both G.-andBinto Eq-. (E-4-9) and taking the Laplace 
transform we obtain . 


( s +2k x s + k x + ,J x^ 8x + ( k xy _ ' > ^ 8y + - (ky Z '“P )#z 

= at, (8a -r(2k + — )6v, + k- 8v. + k 8V. , 

x v x p • dx xy dy xz dz 

(E-51a) 

(k^-vjsx + (s 2 + 2k y s + k^. + p y )6y + (k yz -T)8z 

k* ' 

= J [8a y + V v dx +(2k y + r )8V dy + V 8 v dz ] 


(E-51b) 
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(k zx -p)6x' +(k zy -T)5y +(s^+2k z s + .k^ +.p z )6z 


[6a + k 5v, + k 6v, +(2k + — )6v, ] 

z zx dx zy dy ' z p ' dz 


The' characteristic equation of the system (E-gl) is 


(E-51c) 1 


s +2k s + k’ + u - 
x .x x 


k .— v 
xy 


D(s) =• k yx - v 


+2k s + k’ + u k- - t 
y ; y p x yz 


k zx ‘ P 


k - t 

. zy 


+2k- s . + k 1 + u 
z . z r a 


Equation (E-52) is .identical -to Eq . (E-29) . -In the discussion 

of- Eq. ,(E'29) the conditions that k^j k^> k z , k^ s k^q k^ must satisfy 

were indicated and a method of generating k , k , k was. suggested. 

xy - yz zx 

Then the transfer matrix H takes the- form: 


2k + — 
x p 


k* 

2k ’+ 

y P- 


2k + — 
z P 
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The possibility of a solution which does riot require variable 
-gains was demonstrated in the previous section. It is possible to 
apply the same idea here. 

Let us- choose an H of the following . form: 


H . = 


2k + — 
x p 


,k' 

2k, + -£■ 

y. , p 


k' 

2k '+ — 
z . P 


, (e-54) 


Introducing this H into Eq.' (E-4-9), expand it and taking 
Laplace transforms we obtain. 

[ri 2 + 2k x s + k^ + p x ] 6 x - v6y - p8z -X [8a x +(2k x + F )8V dx- ] 

(E-55a) 

2 -P k' 

- v8x .+ [s- + 2k y s + k^ + p y ]$y t&z '=J^ [6a y +(2k/^£)5 v. ]‘ 

( 

,(E-5?b) 

k* 

- p.8x - r8y +is 4 .+ 2k z s + k| + p z ]8z = £ ^.+(2^+^ )6v az ] 


(E-55c) 


The characteristic equation of 'the system (E -55) is: 


D(s) = 


s + 2k s + k'+p 

X ■ X *x 


's 2 + 2k s + k’+p -T 

y y y 


s + 2k ,s +k'+p , 
z z H zl 


(E-56) 
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. Letting 

2k = Ek = 2k = 2£Sl ana k' = k! = k' = fl 2 (E-57) 

x y . z * j * . ■ 

1 we obtain 

(s £ + ££0s + 0 2 ) 3 + (n x --+ U y + (i K ) (s 2 + 2£SS + 0 2 ) 2 
+ (p- x H y + l* y P z + P x ' .V 2 * P 2 “ .t S ) U 2 + 2 S 0s + A ) ' 

. + P x P y t* z " .^x 1-2 P y P 2 " P z ^ " 2y P T ' = .E-58-). 

Clearly Eg.- (E-58) is identical to Eg.’ (E-42) of the previous _■ 
section; hence the stability question concerning 'the: present .system has 
been- answered there.- 

Hote, however, that the transfer matrix. is different for each 
case and therefore- it operates on 8R. and 6v.. in a different. manner; 
this point may be. of some significance depending on the characteristics 
of 6E r and Sv^. 
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E-5 ELIMINATION OF THE Brp/ER&EHG SYSTEM ERRORS BY MEANS OF 
ALTIMETER INFORMATION 

The' altimeter is a device. which can; measure the. distance of 

the veniele from a celestial hody. 

The barometric altimeter measures the altitude based on 

atmospheric density measurements . 

The radio altimeter operates on the basis of radar principles. 

| . 

The optical altimeter measures the distance of the! vehicle from 

* • " ...... , * . 

a .planet by measuring the visible angular dimension of .the .planet. 

Both tthe radio altimeter . and the optical altimeter' find appli- 
cations in space navigation. The 'barometric altimeter' can .be 'used for 

t . , 

earth bound navigation where the flight takes place inside the atmosphere. 


The idea of using an altimeter for stabilizing the navigation 
system stems from Eq.(D-24) which is: 



(d-24) 

We noted there that the first term in Eq. (d- 24) represents; the' 
gravity errors j.due to tiie position errors included in the numefitor s 
of Eq. (D-3)' and that .the second term of Eq. (D-24) represents the 
gravity errors due to the position errors included in the denominator 
of iEq. (D-3). 



' 95 - 


The altimeter method is. based on 'this observation . The" 

components , s , s , s of s are obtained from the output of the system 
x y z 

and the- altimeter is used- to measure the distance of the . vehicle from 

the center of attraction. If the altimeter measurement includes' an 

error 6h then the output of the altimeter is E + 8h. 

The information E + 6h and s , s , s is fed into the gravity 

x y z 

computer which now simulates the equation 


^XC 


s 

. X 

s yc 



(R+6h) 3 

s y 

g zc 


s 

Lz J 


as contrasted to the equation 


— — 


— , — 

s xc 


s x 

g 

ym 

s 

ye 

s 3 

y 

g 


s 

Ltzcl. 


_zj 


(E-60) 


which is simulated by the gravity computer in all the previous methods 
of damping. Apparently both computers generate identical gravity 
components if. their inputs are error-free; this, of .course, is a basic 
requirement for all gravity- computers . 

The input-output error relations- corresponding to Eq. (E-6 o) are 
Eq. -(D-24) . In the following we develop input-output relations corres- 
ponding to Eq.- (E-59). These error relations will be different frc_ 
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Eg.. (D- 24 ) due to the difference in the.' denominators in Egs. (E- 59 ) 
and (e- 60 ) . Of course we. anticipate" that such a change will affect 
•the stability of the entire system in a favorable manner. ' 

Figure E -6 shows the mechanization of a navigation system 
which uses altimeter information. 

How Eg. (e- 59 ) ean^be' expanded in 'terms .of ■ the errors 6x, 

8y 3 8z and 8h . Consider for example the first of Eg . (E- 59 ) >• this" is : 


SxC = “ ^ (R +8 h )3 


.. x+8x 

V“. 5 

_(R+Sh) 3 


— — - — y = (E+8h) 3 = R 3 jj + higher -order terms in 8h 

(R+Shjf 


(R+8h) 3 


= (x+8x) (R+8h)' 3 = x R "3 + 8 x R -3 - 3 x R""* 1 " 8h + 


higher, order terms 


Introducing Eg. (E-63) into (Eg. (E-6l) we obtain: 


■=27 -8h -+ higher order terms 

R? 


From Appendix D we. have g = - vm and ^ = § = 

x _o „3 R 


This means that Eg. ,(E-6 o) can be rewritten as- follows: 

2 2 3 x 

g xc - = s x " w s 8x + “s R 8h ' + hitler order teittns . 




Figure E-6 Altimeter Stabilization 
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or after dropping Higher order terms 



Sg x 

s xc - 4’ ’■-'•'V 8 * 

. 2 3x 

+ »sF 

6h 

(E-67) 

Similarly 

6g y. " 

gyc - s y = ml 5y 


6h 

(e-68) 

and 

s % = 

g zc - s z = -4 6Z 

+ " 2 3z 
®s E 

8h '• 

(e-69) 


Equations (E-&7)> (E-68) and (E- 69 ) can he written' in. matrix form to 
obtain: - 


<0 

'-2 

-»s 

0 

0 


8X' 


X 

E 

6s y 


0 

2 

-“>s 

0 


8y 

-+3m§8h _ 

1 

S 

Z 

E 

8s z 


.0 

o' 



6z 

z 

- E n 


(E-70) 

Equation (E-70) was ’ anticipated . from 'the disucssian of Eq./ (D 

It is obvious from Eq. (E-70) 'that there - is no coupling in 
error equations for. thq. gravity computer. 

By inspection of Fig. E -6 the - error equations of - the navigation 
loop are: 

6E = sa + 51 (E-71) 

and introducing (E-70). into (E-7l) we obtain: 


8x 

+ ®s 

Sx - = 

• 5a + 
X 

3®s 

X 

R 

8h 

(E-72a) 

6y 

+ ®s 

8y = 

‘ 5a . + 

y 

3»s 

£• 

R* 

8h - 

(E-72b) 

8z 

+ ®s 

8z =, 

8a z ,+ 


£ 

R , 

8h 

(E-72c) 
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Assumi n g that the coefficients o.f Eq. (E-72) vary slowly 


we can take Laplace transforms to obtain 

(s 2 + m 2 )8x = 0 ^[8a x + ’3 fl 2 | 6b] (E-73a) 

(s 2 + to 2 ) .8y = ^[6a y +3o)5 | 8h] (E-73b) 

(s 2 + m 2 ) 8a = <£[Sa z + 3»g | 6h] (E-73c) 

The characteristic equation of the system (E-73) is: 

(s 2 + ro 2 ) = 0 (E-74) , 


Hence the errors vary sinusoidally with, bounded amplitude . - The error 
8h of the altimeter measurement acts as a' forcing. function; its 
effects on the position errors will- be ' small if it can be made small. 

The use of the' altimeter did not stabilize the system 
absolutely; ' it did, however, eliminate the time increasing errors 
and this is' of considerable value'. 
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APPENDIX F 

DIRECTION COSINE EQUATIONS AND THEIR SIMULATION ON 
A DIGITAL COMPUTER 

Alfred R. Schuler 


1- THE DIRECTION COSINE EQUATION 

In order to understand the navigation mechanization thoroughly, 

it is necessary to have a knowledge of the direction cosine equations 

and how they arise. In the 1310011 diagram of Fig. 1, the output 

quantities are A , A , A , SI , IL and a, * which are linear 
V lc v 2c V 3c V lc V 2c T 3 c 
accelerations and angular velocities of the vehicle resolved along the 

vehicular system hut measured with respect to the inertial system. 

The quantities , A^ and A have included in them the gravity 
lc 2c v 3c 

terms . Now in order to determine position with respect to the inertial 
system, it is necessary to know^the orientation of the vehicular system 
with respect to the inertial system at every instant of time. As 
discussed hy Krishnan^, it is necessary to relate points in one coordinate 
system to another rotated arbitrarily with respect to it. Given two 
sets of axes [i] (with components 1^, I £ and 1^) and [v] (with compo- 
nents V^, V 2 and V^), that are arbitrarily oriented it is possible 
to specify the components of one in terms of the other. It is done 
by a series of rotations in a specified order . . Then 


The c's indicate actual .quantities (measured or calculated) at the 
output of the accelerometer loop transfer function. 




IM 

\\ 

Ij = R(cy)B(s)R(Y) 

V 

\V 

W 
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(F-l) 


where E(of)j R(p) and R(y) are rotation matrices. This equation can 
be expanded by multiplying the rotation matrices . The rotations 
specified here are identical to those indicated in Goldstein' 1 ' 




j costyeos0-eos9sin0sinty 

cosi]isin0+cos9cos0sini]i - 

sintysinO ^ 

j\\ 

?I 

= 

-sintycos0-eos9sin0cosi|r -sin4sin0+cos9cosjfcosty 

i 

costysinO 

j 

[ | 

W 


‘ sinQsin^ 

-sinOcosfi 

cose / 

V 


(F-S) 

The elements of the matrix product [d] = R(or)£(g)R('y) represent 
the direction cosines of the angles between the three coordinate axes 
V . , V-g and and the original inertial axes 1^, and 


l h \ 


/ & ii 
1 

d 12 

d 13 ^ 


r\ 

>1 


i 4 21 

d 22 

a 23 


I*') 

W 


\ d„, 
31 

d 32 

d j 




(F-3) 


or [I] 


[D] [V] 
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How it is easy to show that R^y) = R(y)^ 

R-'Ha) = R(a) t (f-4) 

r _1 (b) ' = rCb)* 

That is, the inverse of a rotation matrix is simply equal to its 
transpose - in other words r(q!), R(b) and R(y) are orthogonal matrices. 

Since [d] is the product of three orthogonal matrices, it too 
is orthogonal, that is: 

[D]- 1 = ID]* = Elv^RlplSw* 



(F-5) 


Let Ry be a column matrix whose elements are the components of 
R in the orthogonal coordinate frame V . Components of this same 
vector R in another orthogonal hut rotated frame I are related to by 

*1 = ' D r (F - 6) 

D* is the direction cosine matrix of the V axes referred to the 

V 

X axes. It will normally he used synonmously with D. Also, D^ will 
represent d\ the transpose of D. 

Differentiating Eq. (f- 6) with respect to time yields: 

*1 = D V 


(F-7) 


= BiV + DRy 
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■In general, hoth D -and R are functions of time . 

The time rate of change of a vector R in'-fhe inertial frame is- 
related. to the time rate of change in the vehicular frame hy the 
-.Coriolis operator equation: 1 

/. dR \ _ / dR \ , „ -cr 

( )j “.(■'Jjr \ " “V x w-vj 

How define the matrix [u>y] to he a nine component symmetric matrix 
whose elements are the components r of angular velocity along the axes 
of the V frame and with respect to the I frame. 


■ 0 

-V 3 

“V 

«v 3 

0 

■“V 


-Vi 

0 


w = 


The Coriolis equation can then he written in matrix form: 

Rj = Ry + [my]Ry (F-10). 


Let the unit vecotrs in the two coordinate systems he 

= V jW 

e v = V + "V + \ 

Then 9^ and 9y are the associated column matrices of their components' 

' Then using Eq. (F,-7) 

®I = 0 = + ?% > 

since ,’the unit vectors in the inertial frame are non -rotating constants.’ 
Also using Eq. (F-10). 

'®j. = 0 = 9y + 

Combining these two equations gives 

Dtojytey = My 

D = Dtuoy 3 


or 


(F-ll) 



Expanding the terns in the matrices gives 


d ll 

d 12 

d i3 

l* 

« 


d 21 

d 22 

d 23 

* 

f 


d 3l 

d 32 

d 33_ 


d n 

d 12 

d l3 

d 21 

d 22 

d 23 

_ d 31 

d 32 

d 3?- 


0 





o 


“V, 



Equating components: 

a n = d i2 "V 3 - d i 3 «V 2 
d i2 = a 13 - a i;L ay > 

d 13 = d ll °V 2 " d 12 

a 21 = a 22 “V 3 “ a 23 “Vg 

a 22 = d 23 m V 1 " d 21 "V > 
d 23 = d 21 “Vg " d 22 

d 31 = d 3S “V 3 * d 33 “Vg ' 

d 32 = d 33 °V 1 ' d 31 l 

d 33 = d 31 % " d 32 "V. 

^ 1 / 


(F-12) 


(F-13) 


(P-lA) 
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These equations can he implemented and [D lean he evaluated as 
a function of time. In a digital computing system [uy] is not 
continuously available since the accelerometers are sampled at 
discrete instants of time which are separated by an interval aT. 

Between samplings the [v] frame will rotate and since the angles are 
non-eommutative, computational errors will arise. 

2. SOLUTION OF DIRECTION COSINE EQUATIONS USING DIFFERENCE EQUATIONS 

The digital solution is based upon the replacement of the 
differential equations by difference equations . 

For the first set of differential equations (F-12) we have: 

<*ll(i +1 ) = [a i2 (i) ay (i) - a i3 (i) ny^(i)] &T + d 11 (i) 

d 12 (i + l) = t4 13 (i) ^(i) - c^^i) ^(1)3 fiT + d 12 (i) (F-15) 

a i3 (i + l) = [a iL (l) (Dy^i) - d 12 (i) (i)3 AT + 4^(1) 

The arguments i and i+1 imply samples at times T^ and T^ +1 
respectively. Also, T. + -^ - T. = AT. 

Given 4^(1), d lg (i), d 13 (i)j uy (i)> uy (i) and uy (i), 
these equations can be solved for d 11 (i+l), d lg (i+l), and 
d 13 (i+l) . It is assumed that the inertial package provides the 
necessary angular information at each sampling interval i. 
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Hearranging slightly and putting into matrix form gives: 


a ii (i+1) ~ 


i ay (Da* ^(iJat 


d n(i) 

a i2 (i+1 ) 


-«V (i)at % (i)at 

3 1 


a 12 (i) 

_d 13 (i+l)_ 


(U.. (I)aT. -m. (i)AT 1 

2 l 


- a 13 (i)- 


a xl (i+l) 


1 

P* 

H 

H 

H* 1 


i 

( ± )at 

3 

-«V (1)aT 
2 

d 12 (i+l) 

= [i3 

a 12 (i) 

where L = 

-»V (i)AT 
3 


U)„ (1)aT 
V 1 

d 13 (i+l) __ 


_ a 13 (i)_ 


- «v (i)AT 
2 

-ny (i)at 

1 


The difference equation form for Eq. (F-13) is: 


a 2 i( i+1 ) = [ a 22^ “V &T + a 21^^ 

a 22 ( i+1 ) = [ a 23 (i) (^(i) ‘ a 2 i(i) tOy-^(i) 3 AT + d 2£ (i) (F-17) 

d 23 (i+l) = [d 21 (i) oy^i) - d 22 (i) <y.(i)] AT + d 23 (i) 


Put into matrix form, this becomes: 


a 2l( 1+1J 


\ 1 (i) 

a 22 (i+l) 

= [L] 

d 22 (i) 

d 23 (i +1 )_ 


_ a 23 (i) 
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Likewise* the difference equation for the differential equations 
(F-l4) are 

d 31 (i+l) = [d 32 (i) (1^(3.) ' ^(i) “V A T + ^U) 

d 32 ( i+1 ) = ra 33 (i) 0^ (1) - a 31 (i) (i)] at + a 32 (i) (F-19) 

a 33( :i+1 ) ~ ^31^^ “V " a 32 (i) ujy (i)3 AT + d 33 (i) 



It is to "be noted that [Li is not an orthogonal matrix = This, 
is true even for the case in which rotation is about only one. axis. 

Since the fundamental requirement for a coordinate transformatipn is 
that the transformation matrix he orthogonal, we see that Eqs. (F-l6), 
(F~l8) and (F.-20) will deviate from their desired true values with 
the passage of time . Wo mention has yet been made regarding the 
effect of errors in the measurement of 5^. Errors in this quantity, 
which will be considered as our study advances, will also contribute 
to misorientating the coordinate system- 

The v difference equation approach introduces not only mis orientation 
errors but it also results in a transformed coordinate frame- for which 
the "unit" vectors are n'either mutually ortho'gonal nor of unit length-. 
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In other words , a right handed cartesian coordinate system must satisfy 
the following 'mathematical 1 constraints » 

d i d i = 1 ' d i • = 0 

jj • Jj =' 1 ij - kj. = 0 (F-21) 1 

£ C *5 I = 1 jj * kj = 0 

* V = 1 

■V "V ~ d 

V \ = 1 

ij, jj. and k^. are the unit vectors in the iner.tial frame, 
iy., and k^. are the unit vectors in the vehicular frame, 'Expanding 1 . 
Eq. (F-21) gives: 


d 2 

• 2 2 
+ df^ + = 1 

d., _ &_ + d., ^d + d., „d_~ 

= 0 

11 

12 13 

11 21 12 22 ' 13 23 


d 21 

+d 22 +d 2 3 = 1 

d ll d 3l + d 12 a 32 + d 13 a 33 

= 0 

d 2 i 

+ d 32 + d 33 d 

d 21 d 31 + d o 2 d 3 g + d 23 d 33 

= 0 


*V ' 0 

* \ = 0 ,(F-22) 

= 0 


(F-23) 
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Expanding Eqs. (F-22) yields: 


d ll 

+ d 21 

+ d !l = 

l 

d ll d 12 

+ a 21 d 22 + d 31 d 32 

6 

d L 

+ d 22 

+ d 32 = 

1 

d ll d 13 

+ d 2l d 23 + a 3l d 33 = 

0 

d l 3 ' 

+ d 23 

+ d 33 

1 

d 12 d 13 

+ d 22 d 23 + d 32 d 33. = 

0 


(F-24) 


All of the twelve equations (F-23 and F-2%) must he satisfied 
hy any real orthogonal transformation matrix. 

A usual procedure for evaluation of the direction cosines’ is to 
compute six of them using, for examplf, Eqs . (F-l6) and (f-18) . 

Using Eq. (F-l6), the three components of i-j- are computed in terms of 
iy, jy and k^.. likewise, using Eq. (F-l8), the three components of 
are computed in terms of the vehicular' frame unit vectors. Then, 
the following three equations, taken from Eq. (F-24) are necessary' 
and sufficient to insure that i^ and are othrogonal and each of 
unit length: 


2 2 2 
a U + a 2l + d 3l = 1 


.2 . .2 . ,2 

d i2 d 22 d 32 


(F-25) 


d ll d 12 + d 21 d 22 + d 31 d 32 = 


0 
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The computer must' evaluate the left hand sides' of Eq. (F-25)-. 

If the results do. not- coincide with the desired values, the first 
equation is normalized to unity. 

An algorithm is then provided for satisfying the second two 
equations utilizing the values of d^, d 21 and d^ just computed. 

The three components of TL^ are then evaluated through the relationship 

\ = r i x h 

^ d ll *V + d 12 \ + d 13 x ^ d 21. \ + d 22 + d 23 M 

( d 12 d 23 ‘ d 22 d 13^ \ + ^ d 13 d 21 " d 23 d ll^ 

+ ^ d ll d 22 " d 21 d 12^ \ (F-26) 

3- SOLUTION OF DIRECTION COSINE EQUATIONS BASED UPON A 
TAYLOR SERIES EXPANSION 

2 

Kosmola has suggested a solution in the form of a Taylor 
series . Supposing the elements of to he continuous functions 
of time, the direction cosine matrix 'D is expanded into a Taylor 
series about T^: 

g 3 

d(t^ + at) - d(t^) + d(t ± ) at + d(t^) + *d(t ± ) 


(F-27) 



111 . 


From Eq. (F-ll) 

D = D [ujy ] 

D = D[i»y] + Dtuiy] = Dto ^] 2 + D[tUy] =D ( t'ry ] 2 + [rn^J) 

D = D([roy] 3 + 3[<i>y] t-jjy.] + [ajy]) (f-28) 

etc . 

Thus, we can write Eq. (f- 28) as: 

D ( T i + A*) = D(T.) ^1 + AT (« V (T i )] + Ml 2 {[^(T.)] 2 + [<I v (T i )]} 

+ iyp- rtdJytTj)] 3 + 3 [(By^)] [ffiyCTiH+tBly^i)]-]. 

(F-29) 

Thus, the new transformation matrix DC^ + at) = D(T. +1 ) 
can he derived from D(T^) "by an infinite series of matrix operations . 

We now make a fundamental assumption that makes the problem of 
finding an exact direction cosine matrix solvable . The assumption is 
that ojy remains constant during each sampling interval AT. Under 
these conditions all derivatives of ojy become zero and Eq. (F-29) 
can be written- 



D(\ + AT) = D(T^) { 1 + [u^)] AT + [^(T..)] 2 


]3 ,.(AT ) 3 + 


[(B y( Ti )P J$L. + 

[iBy^i)^ 1 


= D^) e 


(P-30) 
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Let 


= [e v (T 1 )r^ = [ l o v (T i )] dt 


^ + &T 


(f-31) 


and rewrite Eq. (F-l6) as 


D(T ± + AT) = D (T^ ) e 




(f-32 ) 


where [O^I is a skew symmetric matrix sincetufy] is a\skew symmetric matrix. 


tV 


% 0 

3 


r e v 

3 2 

'V 


■% ®V 0 

2 1 


0?-33) 


2 2 2 2 
Let i = 0 + ®v„ + ®V. 


(F-34) 


1 '2 "3 

and note the following two properties of the skew symmetric matrix!©^: 


tV 

= C-D m 

[Qyl 2 je 2 ” 

(F-35) 

2m+l 

= (-D m 

[© v 3 i 2 ” 1 m = 1, 2,-3.. 


[ V 

,..(f-36) 


Using these two recurrence relations in Eq. (F-32) yield: 
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b[t ± + & t] = dOi^) f i + [e v (T .)3 + [e y (T i )] 2 h 

- l, e v ( T i )3 f if" [ VV ]2 i + [ V T i )] 5 T 

+ £ W 32 4 - [e v ( Ti )3 4 


+ ] 


= d(t ± ) 

„2- , k , 6 

+ d(t.) [e v (T.)] { 1 - fr + jj - fr + ”■ 3 


+ d(t.) [e v ( Ti )3 2 ( + ... } 


= H\) 
d(t ± ) to v (a? i )] 


3 5 7 

r e r + ;T AL + ... ‘x 

* " " 3s + ~ jt + j 


D(T.) t 2 t k ,6 ,8 

Te f 1 “ 2T + fT-ST + BT 

i 

D(Tj[e v (i r )]? — 


[e v (T.)] 

= D(T. ) [ 1 + ■ -- " ■■■■■ sin A + 


[e v ( T i> r 


[1 


} 


cos' £}] 


'F-3T) 



13 *, 


How 



0 

" u V 3 &t 

V* 



V' 

0 

lUy^AT 

(F-36) 


'“Vg 61 

%at 

0 



[e v (T..)J 2 


-«V (&T) 2 -t ^ C&T) 2 

(at) 2 

i>V “V (at ) 2 
i 3 


% (at) 2 

1 v 2 

-ay (at) 2 -o^ (at) 2 
3 1 

(IV “V (at) 2 

2 v 3 


2 

ay ay (AT) 

1 v 3 

ay ™ir (at) 

2 v 3 

^( AT ) 2 -^^!) 2 


(F-39) 


where it is understood, that the <u's are evaluated at time t = T 


Also from Eq. (f-3^) 


1/2 


1 - (ay + ay + uy ) AT 

V 1 v 2 v 3 - 


Expanding Eq. (F-37) in terms of its components yields: 


(e-4o) 


'^ll^i+l^ 

^(^i+l^ 

d 13^ T i+l 

d 2l( T i+l^ 

d 22( T i+l^ 

d 23 <Vl 

d 3l( T i+l^ 

d 32^ T i+l^ 

d 33^ T i+l 
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^ 11 ( T i> 

il 2 (T t ) 

d 13 ( V 

/ 

/ 

a 

-d+e 

f+g 

d 21 ( T i) 

d 22 ' ( ' T i) 

d 23^ T l^ 

X 

d+e 

b 

-h+i 

_ a 3l(^ 

d 32 ^ T x^ 

d 33^ 


-f+g 

h+i 

e 


(F-in) 

Where 


1 - XjM 

f = "V® 
, 2 

1 - x 2 M 

g' = ct^M 

1 - XjM 

h ,= ilI 
V 1 

<*v K 

v 3 

i =■ to M 
c 

a, a M 



X 3 

= 2.2 
= «V + “V 
1. 2 

“a' 

= ' % 

"V 

v 2 

X 2 

2 2 
- % +U V 3 

“b 

= “V, 

«v 

v 3' 

*1 

2.2 
- “V + "V 

2 V 3 

“c 

= % 

“V 

3 


(&T) 2 = (&T) g (l - post) = XL - cos l) 

* { % + “V 2 + “vjHi 1 ) 2 “V 


M = (1 - eos i) 



L 6 

„ ... AT sin l sin l 

N — r~ - 

Equation (F-4l) can be solved on a computer to give the 

exact direction cosine matrix at the future time, T. ... . Effort 

l+l 

is. presently being directed towards a) an error analysis for the 
case in which 5^. is not known precisely but has a specified probability 
density function and b) an extension of the above results for the case of 
non zero angular acceleration between samples . 

References 

1. Goldstein, Herbert: "Classical Mechanics," Addison-Wesley 

Publishing Company, 1959* 

3. Kosmola, Albrecht L.: "Feasibility Study of a Gimballess 

Inertial Space Reference," MIl/lL Report R-27^, April i960. 

3. Krishnan, V." "Design and Mathematical Analysis of Gimballess 
Inertial navigation Systems," Chapter 2, a dissertation in 
Electrical Engineering, University of Pennsylvania, Philadelphia, 

Pa., December 1963. 

Ohlberg, Eugene: "Methods of Strapped Down Navigation," 

Nortronics Division of Northrop Corporation, Palos Verdes 
Research Park, California. 

5. Robinson, Alfred C." "On the Use of Quarternions in Simulation 
of Rigid Body Motion," Aeronautical Research Laboratory, 1958. 



A-29217-3 


CONTROL BOUNDARIES 
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